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without entanglement on an ensemble quantum computer 
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Bernstein and Varizani have given the first quantum al- 
gorithm to solve parity problem in which a strong violation 
of the classical imformation theoritic bound comes about. In 
this paper, we refine this algorithm with fewer resource and 
implement a two qubits algorithm in a single query on an 
ensemble quantum computer for the first time. 
PACS: 03.67. Lx 

Associated with the model of quantum computer JEL2|, 
a variety of quantum algorithm has been proposed M 7 1 . 
In the theoretical view, these algorithms have relevance 
to the entanglement phenomena, the peculiar quantum 
property first identified by Erwin Schroedinger [H , which 
is invoked as the mechanism for the speedup of quan- 
tum computing over their classical counterpart J^] . Until 
recently these algorithms were only of theoretical inter- 
est, as it proved extremely difficult to build a quantum 
computer. In the last few years, however, there has been 
substantial progress based on nuclear magnetic resonance 
(NMR) |l(| . Up to now some simple quantum algorithms 
have been realized step by step on NMR quantum com- 
puter, including Deutsch's algorithm ]Tl|-[T7[|, Grover's al- 
gorithm and ordering find's algorithm p3[ . How- 
ever, a sharp criticism has been proposed by Braunstein 
et al. that NMR experiments have not actually realized 
quantum algorithm because at each time step the state 
of the system can be described as a probabilistic ensem- 
ble of unentangled quantum states ^4|. On the other 
hand, some scientists believe that for a specific quantum 
algorithm the power of quantum computer derives from 
quantum superposition and parallelism, other than en- 
tanglement 25 - 28 1 . 

The problem we considered in this paper is the parity 
problem about a database A that contains an arbitrary 
n-bit string a. The answer to queries represented by n-bit 
string x to the database is the parity of the bits common 
to x and a given by (a, x) — a ■ x. Note that the problem 
is to determine a in its entirety, not to merely deter- 
mine the parity of a. the classical determination of n-bit 



string a requires at least n query operations (since n-bit 
string a contains n bits of information and each classical 
evaluation of query operation yields a single bit of in- 
formation). Bernstein and Vazirani have given the first 
quantum algorithm in which n-bit string a can be deter- 
mined in only two queries to the database ||. But by 
preparing the output one bit register in an initial super- 
position (|0) — |1)), the algorithm can be simplified to 
comprise a single query |2^| . Terhal and Smolin rediscov- 
ered this algorithm, which was underappreciated, to solve 
binary problems and coin-weighing problems effectively 
{§■ 

In this paper, based on Bersntein- Vazirani's parity 
problem Q| , we proposed a scheme to solve this problem 
by using less physics resource than the previous algo- 
rithm pj| , [2l| , but without loss of effectivity. Further, we 
demonstrated this algorithm on a two-qubit NMR quan- 
tum computer for the first time. 

Bernstein and Vazirani's parity problem can be de- 
scribed as a function / : {0,1}" — > {0,1}, which is of 
the form f a (x) — a ■ x = (X^i=i a i x i) mod 2, where n-bit 
strings a,x £ {0, 1}™, cii and Xi are the i th bits of a and 
x, and a ■ x denotes the bitwise AND (or mod 2 scalar 
product) a ■ x = (aiAxi) (g> (c^Aa^) <£>■•• <g> (a ra A:c„), the 
answer is to find n-bit string a. The previous quantum 
algorithm [^,^9) to solve this problem theoretically by a 
pair of registers (x, b) , where x £ {0, 1}" , b £ {0, 1}. The 
quantum network to implement the algorithm is shown in 
figure 1, the n + 1 qubits register (x, b) start in the state 
\x) \b) = (|0)) n ^ (|0) - |1)). The function /„ (x) = a- x 
is designed within a unitary operator U / which denotes 
the transform 
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\x) \b) ^ \x) \y®f (x)) = \x)\b®(a-x mod 2)) 
The Hadmard Gate H denote the transform 
|0>-^(|0> + |1» 
|1>->^(|0>-|1» 

If we apply n Hadmard gates = H ® H < 
in parallel to n-qubit, then the n-qubit state transforms 

2™-l 

as ffM : \x) -> -k= Y, (-1)^ \y). Therefore, acting 
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on the input (|0)) n ^ (|0) - |1)), with /„ (x) = a ■ x 



and 



t^t J2 (~ I)" X X y = day, we can evaluate the state of 

y=o 

the input register as 

(|o»*|i> H ^-h IM*)®^ 
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^ E E (-ir x+x - y \y) 



\a) <g> 



x=0 i/=0 
|Q> — 1 1> 



|Q> — 1 1> 
x/2 



It is obviously that one could execute the quantum net- 
work once and measure the n-qubit input register, find- 
ing the n-bit string a in function f a (x) = a ■ x = 
(Et=i a i x i) mod 2 with probability one. 

Our refined version of Bernstein- Vazirani algorithm 
uses n qubits, rather than n + 1 qubits, to find n bits 
string a. A quantum circuit for this refined quantum 
algorithm is shown in figure 2. To compare with the 
original algorithm shown in figure 1, the one-qubit work 
register b is removed because it is redundant in the sense 
that its state does not change. To do so, the binary func- 
tion f a (x) = a ■ x encoded in a n + 1 qubits unitary 
transformation U f was changed into the n qubits propa- 
gator U a such that 



I*) - (-1) A(X) I*) 



(4) 



, and the unitary transformation U a can be decomposed 
as direct products of single qubit operators 



U a = U 1 <g> U 2 



U l = { 



a, 
a-, 
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(6) 

The first Hadmard Gate = H <g> H <g> • • • <g> H 

takes |^o) = (|0)) n to |^> = -4- E I*)- After 

V x=0 

the unitary transformation U a responds to this quan- 

2 n — 1 

turn query, the state is = E ( — 1) 
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The final Hadmard Gate out- 



v x=0 
puts |^ 3 ) 
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\y) = 



Where- 



upon measuring the whole n qubits register identifies a 
with probability 1 (the output states for different a's are 
orthogonal) . 

Quantum entanglement and quantum interfere are usu- 
ally thought to be the key gradient in quantum algorithm 
and the reason why quantum algorithm exceed classical 



algorithm. But in the above refined quantum algorithm, 
there is no entanglement in it. The initial state is \ipo) = 
(|0))", which is obvious separable. After the Hadamard 
transformation, the state is \tpi) = -7= (|0) + |1)) ® 

(|0) + |1)) ® • • • ® (|0) + |1)). Performed by query opera- 
tions U a , the state becomes ^2) = -7= (|0) + e t7Ta ° |1» (g) 

(|0) +e inai |1» O ••• <g> (|0) + e l7m ™ |l)).And the states 
after the second Hadamard transformation is the out- 
put state ^3) = |ao) \a\) ■ ■ ■ \a n ). In the whole proce- 
dure, the state is tensor products of the states of the 
individual qubits, so it is unentanlged. And because 
the operators in the algorithm {H^ n ' , U a and H^) 
are also tensor product of the individual local opera- 
tors on these qubits; = H <g> H ® • • • (g> H and 
U a = U 1 ® U 2 ® • • • (g> U l ® ■ ■ ■ ® f/™" 1 <g> [/". Such a 
unitary transformation cannot change the entanglement 
of a state. 

Experimentally, this quantum algorithm without en- 
tanglement was implemented using the nuclear spins 
of the two hydrogen atoms in a deuterated cytosine 
molecule. |0)(|1)) describes the spin state aligned with 
(against) an externally applied, strong static magnetic 
field Bq in the z direction, the reduced Hamiltonian 
for this two-spin system is to an excellent approximation 
given by 



2tt Jlfl? 



(7) 



where the first two terms describe the free precession of 
spin A and B of two hydrogen atoms about Bq with frc- 
500Mhz, and the chemical shift 



"J A 



2tt 



quencies 2 

t^t I = 765Hz enable us to address each spin (act 
tA 



I 2tt 



as qubit) individually. 1% is the angular momentum oper- 
ator in the +z direction for A, The third term describes a 
scalar spin-spin coupling of the two spins of J ~ 7.17Hz. 
As we know, pulsed radiofrequency (RF) electromagnetic 
fields, oriented in the x — y plane perpendicular to static 
magnetic field B , selectively address either A or B by 
oscillating at frequency u>a and u>b- For example, a RF 
pulse along y rotates a spin about that axis by an angle 9 
proportional to 9 tP, the product of the pulse duration 
t and pulse power P. In this paper, we shall let Ry (9) 
denote 9 rotations act on spin A about y, R AB (9) denote 
9 rotations act on spin A and B about x simultaneously, 
and so forth; superscripts will identify which spin the 
operation acts upon, subscripts denote which axis an RF 
pulse rotates a spin about. 

Experiments are conducted at room temperature and 
pressure on Bruker Avance DMX-500 spectrometer in 
Laboratory of Structure Biology, University of Science 
and Technology of China. A quantum circuit for im- 
plementing this algorithm on a two qubit NMR quan- 
tum computer is shown in figure 2 with n = 2. In 
our experiment, pairs of Hadmard gates were replaced 
by an NMR pseudo-Hadmard gate h (a 90° rotation) 
and its inverse hT 1 pi| . An input peseudopure state 
0o = 1 00) was generated using the approach of Cory 
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et al. HH). This is implemented as R B (tt/3) - G z - 
R% (tt/4) ~t~R^ (-tt/3) - G z , to be read from left to 
right, where G z is the pulsed field gradient along the z 
axis to annihilate all transverse magnetizations, dashes 
are for readability only, and r represents a time interval 
of 1/ (2 J) w 69.735ms. 

The pair of pseudo-Hadmard gates ft and /i" 1 could 
be easily implemented by two hard pulses denoted as 
R AB (tt/2) and R AB (tt/2), the typical pulse lengths 
were 10 — 20us. All unitary transformation U a corre- 
sponding to the query of four possible 2-bit string a — 
{00, 01, 10, 11} in function f a (x) — a-x could be denoted 
as C/ 00 = I A <8> / s , Uoi = / A (8 of, tfio = ^ ® I s , CT U = 
cr^ 4 (8> erf. The C/oo transformation corresponds to the 
unity operation or "do nothing". Uoi and Uw transform 
are separately achieved by applying R B (tt) rotation se- 
lectively on the second qubit B and R A (tt) rotation on 
the first qubit A. The soft z pulse was implemented by 
the time evaluation under the Hamiltonian of Eq. (6) 
with refocusing tt pulses applied at suitable times during 
the evolution period. Since the refocusing tt pulse has the 
effect of time reversal, it can be used to make one term 
in the Hamiltonian evolve while the other terms "freeze" 
Q,^5|. In our experiment, we extended these method to 
realize the soft pulses R A (tt) and R B (tt) separately as 

n/4 - R B (tt) - n/2 - R B X (tt) - n /4, u A Ti = tt 
r 2 /4 - R A (tt) -n/2- R B X (tt) - r 2 /4, uj b t 2 = tt 

the axes of successive tt pulses were chosen in the way 
to cancel imperfections of soft pulses. U\\ corresponds 
to a 7r-rotation R AB (tt) about the axis z of both qubits, 
up to a global phase factor. Global phase changes are 
not detectable in NMR and are hence ignored for the 
purpose of experiment. This non-selectively z- rotation 
was implemented using a composite-pulse sandwich, as a 
set of x and y axes R AB (tt/2) - R AB (tt) - R AB (tt/2) . 

The result of this algorithm are shown in figure 3. Five 
aspects are shown: a reference spectrum acquired using a 
single pulse R y (tt/2), and spectra acquired from the same 
computer implementing the query algorithm for each of 
the four possible functions, f a . Each spectrum consists 
of two closely spaced pair of lines: each pair of lines cor- 
responds to a single qubit, while the barely visible split- 
ting within each pair arises from the spin-spin coupling, 
J. To improve the appearance of the spectra the final de- 
tection pulse was proceeded by a magnetic field gradient 
pulse, which acts to dephase the majority of any error 
terms which might occur. The reference spectrum (a) 
corresponds to the computer being in state |00), and the 
phase of this spectrum was adjusted so that both lines are 
in positive absorption phase(that is, pointing upwards). 
The same phase correction was then applied to the other 
four spectra, allowing positive absorption lines to be in- 
terpreted as qubits in state |0), while negative absorption 
lines can be interpreted as qubits in |1). The left hand 
pair of lines arises from the first spin, and thus corre- 
sponds to the first qubit, while the right hand pair of 



lines corresponds to the second qubits. Thus, spectrum 
(b)-(e) correspond to \a) = |00), \a) = |01), \a) = 1 10) 
and | a) = |11) respectively. It is clear that our implemen- 
tation of a two bit parity problem using a single query 
leaves the computer in a final state, much as expected. 

To summarize, we have presented a refined version of 
Bernstein and Vazirani's quantum algorithm to solve the 
parity problem with fewer resource and implement this 
algorithm in a single query on an ensemble quantum com- 
puter for the first time. This algorithm, which gives 
a strong violation of the classical information theoretic 
bound [ p0| and a clear separation between the quantum 
and classical difficulty of the problem reduces the 
number or queries all the way from n to 1. It is obvi- 
ously that in this algorithm there is no cither entangle 
state or entangle transformation but only the concept of 
coherent superposition is exploited, to prepare "in par- 
allel" an input state which is a superposition of all pos- 
sible classical inputs. Our algrithm and its experimental 
realization demenstrate that the superposition principle 
brings about more effective and concise procedure even 
if the entanglement phenomena do not occur. As we 
all know, some quantum algrithms have relevence with 
entanglement [^,0, but some other not p5|-|28[|, so it is 
meaningful to know the role of entanglement in quantum 
algrithms, i.e., the relationship between the entanglement 
and the complexity of algrithm. 
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upper n line corresponds to n qubits of register X, while 
the lower one line corresponds to one qubit of register b. 

figure 2: The (schematic) refined version of the quan- 
tum circuit shown in figure 1, note that one qubit register 
b is removed by changing the unitary operator. 

figure 3: Experimental spectra from our two-qubit 
NMR quantum computer, spectrum (a) corresponds to 
the initial state 1 00) of register x, and act as a refer- 
ence spectrum; while spectra (b)-(e) were acquired from 
the same computer implementing the refined algorithm 
to solve parity problem in a single query, and determine 
each of the four possible 2-bit string a: (b) a — 00, (c) 
a = 01, (d) a = 10, (e) a = 11. 



Figure caption: 

figure 1: A (schematic) quantum circuit implementing 
Bersntein-Vazirani's parity problem in single query. The 
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